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Abstract 

In the present paper, it is found that dark energy emerges spon- 
taneously from modified gravity. According to cosmological scenario, 
obtained here, the universe inflates for ~ 10~^''sec. in the beginning 
and late universe accelerates after 8.58Gyrs. During the long interme- 
diate period, it decelerates driven by radiation and subsequently by 
matter. Emerged gravitational dark energy mimics quintessence and 
its density falls by 115 orders from its initial value 2.18 x lO^^GeV^ to 
current value 2.19 x 10~^''GeV^. Key Words : Higher-drivative gravity, 
dark energy, inflation and late cosmic acceleration. 
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1. Introduction 

Recent astronomical observations have revolutionized the understand- 
ing of cosmology. Consequences of combined analysis of la Supernovae 
(SNe la) observations [1, 2], galaxy cluster measurements [3] and cos- 
mic microwave background (CMB) data [4] has shown that dark en- 
ergy dominates the present universe causing cosmic acceleration. This 
acceleration is observed at a very small red-shift showing that it is a 
recent phenomenon in the late universe. Observations of 16 type la Su- 
pernovae made by Hubble Space Telescope (HST)[5], modified earlier 
astronomical results and provided conclusive evidence for deceleration 
prior to cosmic acceleration caused by dark energy in the recent past. 

The acceleration is realized with negative pressure and positive en- 
ergy density that violate strong energy condition (SEC) [6]. This vi- 
olation of SEC gives reverse gravitational effect. Due to this effect, 
universe gets a jerk and transition from deceleration to acceleration 
takes place. Of late, phantom energy has appeared as a potential dark 
energy candidate in this arena, which violates weak as well as strong 
energy condition [7, 8]. 

Although dark energy dominates today , its present density is very 
low ~ 2.18 X 10~^^GeV^. The dark energy density, in the early universe, 
is found to have been very high. For example, the density at the Planck 
scale is found to be 1.18 x lO^^GeV^, at the GUT (grand unified theory) 
scale it is ~ lO^^GeV^ and at QCD scale it is ~ IQ-^GeV^. 

The simplest dark energy candidate is the cosmological constant A, 
but it needs to be extremely fine-tuned to satisfy the current value of 
the dark energy density. Alternatively, to explain decay of the density, 
many dynamical models have been suggested, where A varies slowly 
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with cosmic time (t) [9, 10]. In addition to models with dynamical 
A(t), many hydrodynamic models with or without dissipative pressure 
have been proposed in which barotropic fluid is the source of dark 
energy [10]. Chaplygin gas as well as generalized Chaplygin gas have 
also been considered as possible dark energy sources due to negative 
pressure of these [8, 11, 12, 13]. 

In field-theoretic models, scalar field sources of dark energy are the 
most natural. Apart from quintessence and phantom [6, 7, 8], other 
scalar field models involve k-essence [14] and tachyon fields[15]. Ref. 
[16] proposes dark energy models derived from the action with string 
gravity corrections as well as dilaton and modulus fields. 

Other than these approaches, some authors have considered modi- 
fied gravitational action by adding a function f{R) {R being the Ricci 
scalar curvature) to Einstein-Hilbert lagrangian, where f{R) provides 
a gravitational alternative for dark energy causing late-time acceler- 
ation of the universe. In [17, 18], f{R) = l/R, but this has the 
problem of instabilities of gravitationally bound objects [19]. In [20], 
f{R) has the form aR^ + hR~'^ (where a,h,m,n are real numbers). 
The R term with a positive power is dominant in the early universe 
giving power- law inflation for 1 < m < 2. The R term with a neg- 
ative power dominates in the late universe, when curvature is small 
and gives late-time acceleration. In other attempts adopting this ap- 
proach, f{R) has been taken as InR [21] and Tr{l/R) [22]. Some 
other modifications are suggested in [23, 24]. Recently, a function of 
Gauss-Bonnet term has been considered for f{R) as gravitational dark 
energy [25]. In [26], f{R) is taken as Yukawa-like term of R and in [27], 
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f{R) oc {aR'^ + bRt'''Rf,^ + cRf"'P''R^i,i^„y^, which leads to current accel- 
eration without resorting to dark energy. Thus, in the past few years, 
many models have been proposed to explain accelerated universe. A 
review on modified gravity as an alternative to dark energy is available 
in [28]. A more comprehensive review is provided in [31]. 

All of these models are phenomenological in the sense that an idea of 
dark energy is introduced a priori either in terms of gravitational field 
or non-gravitational field. In spite of these attempts, cosmic accelera- 
tion is still a challenge. So, a model is required that can explain early 
inflation followed by a decelerated phase and then late acceleration be- 
ginning in the recent past. Here, an attempt is made in this direction 
explaining fall of DED from a very high value to current low value. The 
gravitational origin of dark energy is discussed using an approach dif- 
ferent from the approach in [17, 18, 19, 20, 21, 22, 23, 24, 25, 28], where 
non-linear term of curvature is considered as dark energy lagrangian. 
In recent work [30], it has been shown that f{R)- dark energy mod- 
els do not give viable cosmology. The model considered here is not 
f{R)- dark energy model, but it is a cosmological model based on 
/(i?) -gravity. There is a crucial difference between these two. In the 
former case, which is criticized in [30], non-linear terms of curvature 
are treated as dark energy. By contrast in the latter case, neither linear 
nor non-linear term is considered as dark energy.In the present model, 
dark energy emerges from linear term (Einstein-Hilbert term) and non- 
linear term aR^'^~'^^ {a is a coupling constant) of the scalar curvature R 
in the action. DE, obtained here, mimics quintessence. Although it is 
induced by the gravitational sector, a corresponding scalar is derived 
here. In the following section, it will be shown that r = 0.75(1 -|- Wde) 
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with Wdc being the equation of state parameter for dark energy In the 
case of quintessence DE, — 1/3 > Wde >— 1, so0<r<l/2 yielding 
2 > (2 - r) > 3/2. As (2 - r) > 0, the term aR^'^~''^ does not cause 
instabihty problem mentioned previously [17, 18, 19]. 

In the present model, the Priedmann equation is derived with dark 
energy , dark radiation and dark matter density emerging from the 
gravitational sector as the combined effect of curvature and a scalar 
field X- Spontaneous emergence of dark energy from the gravitational 
sector, which is a physical concept, is caused by dual roles of as a 
physical as well as a geometrical field [32, 33, 34, 35] . 

This leads to a very interesting cosmology. Here, cosmological pic- 
ture begins with power-law inflation, driven by curvature-induced dark 
energy, for a fraction of a second and universe expands 10^^ times. 
It is shown that elementary particles are created due to decay of the 
curvature-induced dark energy scalar during the inflationary period 
being dubbed phase I. Radiation is emitted due to decay of fluctua- 
tion of at the end of this phase. This radiation heats the universe 
and existing elementary particles attain thermal equilibrium with emit- 
ted radiation quickly. Temperature of this radiation is evaluated to 
be 1.4 X lO^^GeV. This event heralds the big-bang scenario giving a 
possible geometrical origin of cosmic microwave background radiation. 
Thus, we have two kinds of radiation (i) dark radiation and (ii) emit- 
ted radiation due to decay of fluctuations of 0. The emitted also ther- 
malizes dark radiation and radiation component (including elementary 
particles and dark radiation in thermal equilibrium with emitted radi- 
ation) in Friedmann equation attains the temperature 1.4 x 10^''GeV. 
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During the inflationary period, dark energy density falls by 3 orders. 
As a result, radiation dominates and universe decelerates. 

It is mentioned above that dark matter emerges as a combined ef- 
fect of curvature and another scalar x- The dark matter, so obtained, 
is non-baryonic. Apart from it, it is proposed that baryonic matter 
is formed from created elementary particles during phase 1 (inflation) 
through various proceesses of standard cosmology such as nucleosyn- 
thesis, baryosynthesis and hydrogen recombination, which are not dis- 
cussed here. 

According to Wilkinson Microwave Anisotropy Probe (WMAP) [29], 
around 386 kyr of age, matter begins to dominate over radiation and 
deceleration continued . At red-shift Z^^, ~ 0.4, a transition takes place. 
Dark energy re-dominated and decelerating universe got speeded up. 
Dark energy density, obtained here, falls from 2.58 x lO^^GeV^ to its 
current value falling by 115 orders. Thus, this model explains early 
and late cosmic acceleration including deceleration for a long period 
between these two phases consistent with standard cosmology results. 
Thus, this model gives a unified picture of cosmological evolution. In- 
terestingly, it is found below that this is the situation of our observable 
universe, which resides in a larger unobsrvable universe. 

The paper is organized as follows. Section 2.1 contains action of the 
theory, basic field equations and manifestation of curvature scalar as a 
physical field. Friedmann equation, giving cosmic dynamics is derived 
in section 2.2. Scalar field models for dark energy have been very 
popular. Though dark energy has gravitational origin here, it is still 
interesting to probe scalar field 0and its corresponding potential that 
give rise to it. This probe is done in section 2.3. In section 3, Friedmann 
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equation is solved for when dark energy drives the cosmic dynamics 
and it is found that universe inflates with power-law expansion for 
an extremely short period. The creation of elementary particles and 
radiation, due to decay of (f) is also discussed in this section. Moreover, 
this section contains discussion on other three phases of the universe. In 
section 4, idea of larger universe is introduced which is inspired by fall 
of DED by 115 orders. Results of the paper are summarized in the last 
section. 

Natural units {h — c — 1) are used here with GeV as the fundamental 
unit, where h and c have their usual meaning. In this unit, it is found 

that lGeV~^ = 6.58 x IQ-^Sgec. 



2. Dark energy Caused by dual role 
of the Ricci scalar curvature 
2.1 Dual role of the curvature scalar 
Here, investigations begin with the action 

S^J d'x^[-:^ + aR^'-^^ + [^g^-^d,xd.X-V{x)}], (2.1) 

where R is the Ricci scalar curvature, G — Mp'^{Mp — 10^^ GeV is the 
Planck mass), a is a coupling constant with dimension (mass)^^ with 
r being a real number and x is a scalar field with potential V{x)- It 
is shown below that (1 — r) > 0. So, the instability problem does not 
arise. Here, a scalar field x with potential V{x) is introduced, which 
yields dark matter density. 

The action (2.1) yields gravitational field equations 

- Iqi^uR] + a[(2 - r){V^V,i?(^-'^) - 2g^,nR('-^^+ 
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(2.2) 

using the condition SSg/Sg^" = 0. Here, denotes covariant derivative 
and the operator □ is given as 

19/ , — „,, d 



-g dxf^ \ dx^ ■ 

with n^v — 0, 1, 2, 3 and g^t,^ as metric tensor components. 

Taking trace of eqs.(2.2), it is obtained that 

r 1 r 
UR = v^i?V^i? + 



R ^ 3(2 -r)(l -r) lienGa 

-^{-d^xdaX + '^V{x)}-rR' 



(2.3) 



2a 

with q; > to avoid the ghost problem. 

Eq.(2.3) is typically of the Klein - Gordon form with the source term 
given on the right hand side of this equation. It shows that the Ricci 
scalar R behaves as a physical field in addition to its usual role as a 
geometrical field [33, 34, 35, 36]. 

Moreover, from the action(2.1), equation for x is obtained as 

Ox + V\x)-0 (2.4a) 
with V'{x) — dV/dx and components of energy- momentum tensor as 
= d'^xd'x - 9'" {Id'^xdaX - V{x)) ■ (2-46) 

2.2 Emergence of dark energy 

Experimental evidences support spatially homogeneous fiat model of 
the universe [37]. So, the line-element, giving geometry of the universe, 
is taken as 

dS"" = df - a\t) [dx"" + dy' + dz""] (2.5) 
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with a(t) as the scale factor. 

In the space-time, given by eq.(2.5), eqs.(2.3) and (2.4a) are obtained 

as 

■■ rii^ 1 ri?(i+'') 

R + 3-R = — — + 

a R 3(2-r)(l-r) Ll67rGa 

2^ 

and 



.{^x' + ^V{x)}-rR'] (2.6) 



X + 3-X + V'{x)-0 (2.7) 

(X 

as x^^iUi^i't) = xif) due to spatial homogeneity. Here, overdot gives 
derivative with respect to time t. 

In most situations ( matter-dominated model and dark energy dom- 
inated models), we have a{t) as a power-law solution, for which R is 
obtained as power-law function of a{t). So, it is reasonable to take 

= ^, (2.8) 

where n 7^ is a real number and A is a constant with dimension 
(mass)^. 

Using eq.(2.8), eq.(2.6) yields 

d / a 



-y +[3-n(l-r)](-j = 3,(2 -,)(!-,) 



■(^x' + nx))+3^x' = 0. (2.10) 



Eq.(2.7) integrates to 

d (I . 
di 

Moreover, in the space-time (2.5), energy density and p^ are ob- 
tained from eq.(2.4b) as 

Px^lx' + Vix) (2.11a) 
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and 
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(2.116) 



For a cosmic fluid, pressure is proportional to its energy density. So, 
for a cosmic fluid, given by the scalar fleld x, 

Px = kpx {2. lie) 

with k being the proportionality constant. From eqs.(2.11a), (2.11b) 
and (2.11c), it is obtained that 



v{x) 



[1 - k] ^. 



X 



2[l + k] 

Now using k — n{l — r)/3 in eq.(2.11d) 
V(x) = 

2[3 + n(l-r)]^ ■ 
Connecting eqs.(2.12) and (2.10) and integrating, we obtain 



(2.11d) 



X 



B 



a^+n{l-r) 

with B being an integration constant with dimension (mass)^. 
Further, (2.12) and (2.13) yield 



(2.12) 



(2.13) 



-X^ + W{x)^ 



3[1 -n(l -r)] 



B 



(2.14) 



[3 + n{l - r)] a3+"(i-0 ' 
Connecting eqs.(2.9), (2.11a,b), (2.12), (2.13) and (2.14), we have 



3n(2-r)(l-r) 



l^T^GaA^ 



{l-n(l-r)} BA^""-^^ A 
2q;[3 + n(l - r)] ^ ^a"J 

In the space-time (2.5), the scalar curvature is calculated as 



= 6 



dt\a/ \a/ 



R, given by eqs.(2.8) and (2.16a) yield 

d /d\ _ /d\ ^1 A 



6 



dt 



(5) +2(5)1=^ 



(2.15) 



(2.16a) 



(2.166) 
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Elimination of j from eqs.(2.15) and (2.16) leads to the Fried- 
mann equation 

1 



0' 



3n(2-r)(l-r)[l-n(l-r)] 



{l-n(l-r)} BA^'^-^h 1 / r 1\ ^ 

2a[3 + n{l-r)] J ^ [1 - n(l - r)] V3n(2 - r)(l - r) ~ 6/ 

(2.17) 

giving dynamics of the universe. 

Setting n — A, eq.(2.17) is re-written as 

1 



12(2-r)(l -r)[-3 + 4r] 



leTrGaA"- 



{-3 + 4r)}SA('-i)] 



1 / r \A 

2Q;[7-4r] J ^ 6[-3 + 4r] 12(2 - r)(l - r) ~ V ^ ^^'"^^^ 

It is interesting to see that, in eq.(2.18), the second term on the 

right hand side ( being ~ ) has the form of matter density pm ( for 

prcssurclcss baryonic and non-baryonic matter), which is the combined 

effect of curvature and scalar field x with potential given by eq.(2.12). 

As this term emerges due to x scalar in the action (2.1), it is dark 

matter density pdm- The term ~ has the form of radiation density 

emerging from curvature. Such a term appears in brane-gravity induced 

Friedmann equation being termed as dark radiation [38]. So, in an 

analogous manner, here also, the term proportional to ~ in (2.18) 

is recognized as dark radiation density pdrd- The first term, on the right 

hand side of eq.(2.18) being proportional to a~^'^ and emerging from 

the gravitational sector, is recognized as dark energy density pdc- 

It is reasonable to remark here that, in this approach, dark energy 

density emerges due to presence of non-linear term of R in the action 

(2.1). It is unlike the approach in [17, 18, 19, 20, 21, 22, 23, 24, 25, 28], 

where non-linear term of R is used as lagrangian density for dark energy 
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in the action. Here, tlie non-linear term i?'-^^'-' is not dark energy 
lagrangian. 

Now, eq.(2.18) looks like 

/a\ 2 SttG . , X 

\a) ^ "^'^'^de + Pdrd + Pdm), (2.19a) 

where 

SnG ar^ 

" ~ 192(2 - r) (1 - r) [-3 + 4r]nGaA^ ' ^ ^ ^ 

SttG _ 1 / r \A 

^^drd - g[_3 ^ l2(2-r)(l-r) " V ^ ^ ^^''^ 

and 

8nG _ (-l + 4r) 

-^P^"^ - ~24a(2-r)(l-r)(7-4r) " ^ ^^^"^^ 

Conservation equation for dark energy is given as 

Pde + 3-(l+Wde)pde = 0, (2.20) 

a 

where Wde = Pde/Pde with Pde(Pde) being pressure (density) for dark 
energy. 

Connecting eqs.(2.19b) and (2.20), it is obtained that 

4r = 3(l + Wde). (2.21) 

Using eqs.(2.21) in eqs.(2.19b,c,d), we obtain that 
SttG _ _ a-3(i+wde) 

3 36Wde(5 — 3Wde)(l — 3Wde)7r(jQ;A'' ' 

. ^ ^ - l) 4 (2.22.) 



3 18wdeV4(5-3wde)(l-3wde) Ja^ 

and 

SttG^^ ^ 8(2 + 3wde) BA(^+^-'^^^ 

3 ™ 3a(5 - 3wdc)(4 - 3wde)(l - 3wde) 
Eqs.(22a,b,c) show that for positive Pde^Pdrd a-nd pdm 



Wde < -2/3. 
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According to WMAP [29, 39], current values of pdm and pdc are p^^^ = 
0.23p°, and p[J, = 0.73p°, with pO, = SH^/SnG, Hq = lOOhkm/Mpcsec = 
2.33 X 0.68 X 10~^^GeV using h = 0.68 ( a value having the maximum 
hkehhood ). 

These observational values are used to determine constants in eqs. (2.22a) 
and (2.22c). As a result, 

P.. = 0.73p2,(^2)'"""'-' (2.23a) 

and 

Pd„, = 0.23p°,(^)', (2.236) 

where Oq is the present scale factor. 

Now, using eqs.(23a,b), Priedmann equation (2.19a) takes the form 

(3^--o[-K?)"^'"--(?)V¥-- p--) 

First two terms, on the right hand side of eq.(2.24) arising due to 
curvature scalar R, are physical concepts, whereas term on the left hand 
side is geometrical and it is caused by OR in eq.(2.3). All these terms 
arise from the gravitational action (2.1). It is due to manifestation of 

as a physical field as well as geometry[32] . 
2.3 Scalar field for dark energy 

Here, dark energy is obtained from modified gravity. But, in non- 
gravitational models, it is caused by some exotic matter given by a 
scalar field (quintessence, tachyon and phantom). Though, we have a 
gravitational origin of DE here, it is natural to probe a scalar (p (which 
is not a field for an exotic matter) giving p^e obtained above. With 
V{(f)) as potential, (p obeys the equation 

+ 3-0 + y'(0) =0, (2.25) 
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and has mass dimension equal to 1. Here, = dV/dcp. DE density 

and pressure, in terms of 4>, are given as 

Pde = ^0' + n0) (2.26a) 



and 



Pde = - V(ct>) (2.266) 



It is interesting to sec that conservation equation (2.20) yields (2.25) 
on using (2.26a) and (2.26b). 
Eqs. (2.26a) and (2.26b) yield 

Pde = (2.27a) 

1 +Wde 



and 



^(0)^l(l_Wde)pde (2.276) 



using pde = WdePde. Now, from eqs. (2. 21a) and (2.27a) 



and 



n0)^|i(l-w.)i/.^(^)*^"=' (2.2m) 



3. Cosmic evolution 

In what follows, cosmic evolution passes through four different phases. 

Phase 1. Power - law inflation 

According to the present cosmological picture, in the beginning, the 
universe contains DE only. So, Priedmann equation is given as 

= + = (3.1) 

using eqs. (2.26a) and (2.27a). Here H — a/a. 
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(2.24) looks like 
Eq.(3.3) yields 

■ /3(1 + Wde) „ . 

<f>-H^^-^ (3.2a) 



and 

3HH' 1 



inG 2 d0 
using eq.(2.25). Here, H' = dH/dcf). 

Integrating eq.(3.2a), we obtain 



+ V'{(j)) = -3H^ {3.2b) 



a = a(0)e*/ ^^i+^da) (3.3) 



taking 0(0) = without any harm to physics. 
Connecting eqs. (2.28b) and (3.3) 



2 1Q 

m = I^(l-^de)i^o' 



ao 



(3.4) 



l3(l+Wde) _ /4iG 

.a(0)J 

Thus, we obtain dependence of the scale factor on (p during DE 
dominance in the early universe. Using DE dominance, in the other 
form Priedmann equation, given by (2.24), we obtain 



which is integrated to 



2/3(l+Wde) 



(3.6) 

where a(0) = a{t = 0) 7^ 0, which is justified in what follows. 

At the beginning of the universe ,i.e. at t = 0, energy density p^e 
is expected to be very high, but not infinite. According to eq. (2.22a), 
Pde is infinite, if a(0) = 0. It suggests that a(0) should be non-zero to 
avoid the unphysical situation of divergent pde- 
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As Wde < —1/3, 3(1+Wdc) > 2. Hence a{t), given by eq.(3.2), exhibits 
power-law inflation as d{t) > 0. 

Further, eqs.(3.2a) and (3.2b) imply 



H' _ _^ /3(l + w,.) p^^j 



inG V SttG 

giving the slow-roll parameter 

Experimental results restrict e < 0.05 [43]. So, taking e ~ 0.05, 
eq.(3.8) yields 

3(1 + Wde) 1 



and 



20 P-^") 



29 

Wde = (3.96) 



It shows that DE, obtained here, mimics quintessence as Wde > — 1. 
The maximum number of e-fold for the inflation [43] is 

^max ^ g3_3 ^ Q 25ln€ ~ 62.55. (3.10) 
If Oc is the scale factor at the end of inflation 

4£- = 6^2.55 ^ 1 46 X 1027. (3.11) 

a(0) ^ ^ 

According to observational data, given above (in sub-section 2(b)), 
the present DE density = 2.18 x lO^^^GeV^ 

In the present time,pde acquired the value p^^ (given above) falling 
by 115 orders from the beginning , 

p(0) = lO^^Vde = 2.18 X lO^^^GeV^ (3.12a) 

(2.23a) yields 

- n -1 3(1+Wde) 

Pde(O) = -7^ ■ (3-126) 
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Using cqs.(3.9a), (3.11) and (3.12a), (3.12b) yields 

p^^(Q = 4.19 X lO^^GeV^ (3.12c) 
(2.27b) and (3.12c), it is obtained that 

Vc = 4.12 X lO^^GeV^ (3.12d) 

It is interesting to see that Vc satisfies the observational constraint 

[43] 

-2.6x10-9 



for perturbation amplitude at the end of infiation. where K is the 
potential at the end of infiation. Prom eqs.(2.23a) and (3.12a), it is 
obtained that 

ao 1^(^^--)^io-^GeV^ (3.13a) 



.a(0) 

yielding an exceptionally large value of ao/a{0) as 



-10^15° (3.136) 



a(0) 

For the time being, we postpone a discussion on this result to the next 
section. 

Using the results (3.9a) and (3.13a), eq.(3.6) is obtained as 

a{t) = a(0)[l + 1.52 x 10^^^]'°. (3.14) 

Incorporation of the result (3.11) for a* in eq.(3.14), the inflation period 
tc is evaluated as 

tc ~ 1.44 X 10-^^GeV-^ = 9.44 x lO'^^^sec. (3.15) 

At the end of inflation, it is obtained that 



connecting eqs.(3.3),(3.9b), (3.11) and (3.12d) . Here 0c — 0(^c)- 
3.1.1 Pcirticle creation during phase 1 and thermalization 
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It is obtained above that 0, given by (3.16), is the scalar for the 
curvature induced quintessence dark energy. At the end of inflation, 
when = (/)^,\/ = 14 = 4.12 x lO^^GeV^ being given by (3.12d). 
Connecting (2.27b), (3.9b) and (3.12a), it is obtained that at i = 0, F = 
Vq is given as 

Vo = ^{l- Wde)Pde(O) = 2.14 X lO^^GeV^ (3.17) 

These results show that V{(f)) falls by 4 orders, when (p rolls down 
from the state (f) — to (p — (pc- During this period, early universe 
inflates. 

If $(t,x) and V^(t,x) are spinless field and spin- 1/2 Dirac spinor 
present during this phase, it is natural for these quantum fields to 
interact with curvature induced quintessence scalar with interaction 
terms — and — /i'00'0('0 — V'^T^ with being adjoint of ip and 

7*^ being time-component of Dirac matrix in fiat space-time. Here g 
and h are dimensionless coupling constants. 

(i) Creation of spinless bosons 

In the homogeneous space-time (2,5), the equation (Al) (Appendix 
A) for (j) with mass m$ and interaction term with is obtained as 



.a 



-2 



Q2 Q2 Q2 



$(t,x) + [^0'(t)+rr4]$(t) =0. (3.1^ 



- dx^ dip' dz^ - 

Using decomposition (A2) given in Appendix A, (3.18) is obtained 

as 



a • 

$ + 3-$ + 
a 



— + g(j>\t) + mi *(t, x) = 0. (3.19) 



Using the approximation 
a{t) a(0)[l + 1.52 x 10^^^]^° ~ a(0)[l + 3.04 x IQ^H] (3.20) 
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for the scale factor during phase I, conformal time r] is defined as 
rj = J dt'/a{t') = [3.04 x 10^^a(o)]-^Zn[l + 3.04 x IQ^H] (3.21) 
Connecting (3.20) and (3.21), 

a{fi) ~ 0(0)6^ (3.22a) 

where 

f} = 3.04 X 10^^a(o)r;. (3.226) 

Using the conformal time 77 and defining ^%{fi) — a{fj)^k{t), (3.18) 
reduces to 

"^'^^ + [P - 1 + + m|}a2(0)e2^]$^ = 0, (3.23) 



where 



and 



k = A;/3.04 x lO^^a(O), k = A;/[3.04 x lO^^a(O)]^ 



m$ = m$/3.04 x lO^^a(O). 

Now two cases arise (i) when (f) is very small i.e. near to ~ and 
(ii) when (p is large i.e. near to ~ 4>c- 
(a) Case 1.0 is near to ~ 

In this case, (3.23) reduces to 

^ + [P + rhla\0)e^'']^ = 0, (3.24) 

ignoring 1 compared to other terms. 
(3.24) yields the normalized solution as 

T(i + ifc)r(i -ifc)ii/2 



Jf,(m^aiO)e^y (3.25) 



2k{2nf 

The normalization is done using the scalar product defined in (A6). 
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So, it is obtained that 



T(i + i^)r(i -iAi)ii/2 



2A;(27r)2 



(3.26) 



in— and owt— states are obtained when r] tend to —oo and cxd re- 
spectively. Thus 

T(i + ifc)r(i-ifc)iV2 



and 



2A;(27r)3 
r(l + i^)r(l -ifc)ii/2 



«"^(-^)<^ife(^f«(0)e~*') (3.27a) 



a"'(^)^ife(^*«(0)e^). (3.276) 



2A;(27r)3 

Using definition of and from (A8,9) and (3.27a,b), it is obtained 
that 



1 + A;2 



and 



[a2(0)A;]-2 

Condition (AlO) and results (3.29a,b) yield 

a\0)k = 1. 



[am 



-4 



(3.29a, 6) 



(3.30) 



Non-zero value of A^^ = shows creation of particle - antiparticle 
spinless bosons [40, 41]. 
(b) Case 2. is near to p c:^ p,, 

Proceeding in the same way, it is found that, in this case also, results 
for Ictfep and |/5fcp arc obtained as given in (3.29a,b). 

The decay rate of 4>, due to creation of scalar particle-antiparticle 
pairs [40, 41], is given by 

r,^#* = M n l«.r^)M = ^4TX^?^[i + i.52xio^^r'^ 

(3.31a) 
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using (3.29a), (3.30) and 4- volume 

The decay can continue till decay rate T^^^q, is larger than the 
expansion rate H of the universe. So, 



> H. 



It yields 



9.27 X lO^^a(O) , , 3.04 x 10^^ 

^ ^ [1 + 1.52 X IQ^H]-^^ > J- -— — — (3.32a) 



(l + a4(0)) ^ ■ ^ - [1 + 1.52 X 1014^]' 

(3.32a) shows decay of to scalar particles upto the time 



t ~ 6.7 X 10 



-15 



a(0) 



1/60 



-1 + 0^(0) J 

(ii) Creation of spin-1/2 Dirac fermions 



(3.326) 



The action of Dirac field iJj with interaction term —hijj(j)ijj and 
yields the equation 

iYD/^ip - [h(j) + m^]i/j = (3.33) 

where Dn is defined in (All). 

Using the conformal time rj (given by (3.21)), (3.33) is obtained as 

[fdo + fdi + fd2 + fds + i{H + m^)a{r])]^ = 0, (3.34) 

where ^ — a'^^'^{T])ip. 

Here also, wc take two cases as above, 
(a) Case 1. is near to ~ 

In this case, (3.34) reduces to 

[fdo + fdi + fd2 + fd3 + im^a{7])]^ = 0. (3.35) 

Operating (3.35) by [7° 9o+7^9i+7^ 92+7^^3— im^a(?7)] from left and 
using decompositions of given in Appendix A, differential equation 



(3.36a) 



(3.366) 
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for fk,s{9k,s) is obtained as 

/" + [k^ + im^ea'iri) + mla^r])]/ = 0, 

where / = fk,s{gk,s),l^Us = eUs, ^^Us = (^Us and e = ±1. 
Introducing (3.22a,b), (3.36a) is re- written as 

^ + [K' + zm^ea(0)e'^ - mla^{0)e^^]f = 0, 

where K = k/3M x 10^^ and fh^ = -im^/3M x 10^^ 
So, (3.36b) is approximated as 

ignoring im.0ea(O)e'' = (m^/3.04xl0^^)ea(0)e'' compared to (m^/3.04x 
The normahzed solution of (3.37) is obtained as 



(3.37) 



/ 



r(i + /A')r(i - iK) 



1/2 



Jii^(m^a(0)eV3.04 x 10^^). (3.38) 



Connecting (A14) and (A15) to (3.38), it is obtained that 



T{l + iK)T{l-iK)iy^ 



-3/2 3rj/2 



Oq e 



2ir(27r)3 

JiK{m^a{0)ey3M x 10^^)e-'^-^Us, 



X 



(3.39a) 



2K{2Tif 



l'6\Ak.x 



Jiif(m^a(0)eV3.04 x 10^^)6 



(3.39b) 
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Using in- and out-statcs defined as above we obtain ipfj^ ^ and 
for ?7 < and f] > respectively. Similarly, _g-^ and ipjY{-k -s) 

are obtained. 

So, definitions of ak,s and Pk,s, given by (A 17) and (A 18) subject 
to the condition (A 16), it is obtained that 

= I/^mI' = ^ (3.40a, 6) 

which shows creation of spin- 1/2 particle-antiparticle pairs, 
(b) Case 2. cp is near to ~ 0c 

Like spinless case, here also, |q!a;,sP and \Pk,s\'^ are obtained as given 
in (3.40a,b) for cf) c:i (f)c . 

The decay rate of ■0, due to creation of spin- 1/2 particle-antiparticle 

pairs [40, 41], is given by 

v^, . TT , ,2/.. 1.854 X 10^5/n2, 
r^_.v;v = J2 ^< n = [1 + 1.52 X IQ'Hf^ 

s=±l K=-oo " 

(3.41) 

using V4 given by (3.31b). 

The decay of ^ -0-0 can continue till decay rate T^_^^^ is larger 
than the expansion rate H (obtained from (3.14)) of the universe. So, 

It yields that this decay will continue till 

t ~ 6.4 X 10-^^GeV-\ (3.42) 

3.1.2. Fluctuation of at the end of inflation and thermalization 



Here quintessence scalar is playing the role of infiaton causing power- 
law infiation, so hke infiaton, fiuctuates as = 0c(i) -|- 50(i, x) at the 
end of infiation. 
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From (2.25), Klein-Gordon equation for x) is obtained as 

Q2 Q2 Q1 



50(i,x) + 3-(50(i, x) — a 
a 



+ 



(50(t,x) 



+l^"(0e)50(t,x) = O, 
Using 50(i,x) = (f)f{t)e~''^-'^ in (3.43 ), it is obtained that 



(j)f + 3-0/ + 



+ 3.68 X 10 



28 



(3.43) 



(3.44a) 



Using approximated a{t) from (3.20) and r = 1 + 3.04 x lO^H, (3.44a 
)is obtained as 



dV/ ^ 3 #/ ^ 



i;2 

+ 3.98 X 10-^ 



(3.446) 



where /cc = A;/3.04 x lO^^a(O). This equation yields a solution 

(Pf = Er-Vo.998(0.063r). (3.44c) 
Using the approximation 



Jn{x) ~ ^y2/^lxcos[x — TT / 4 + n7r/2] 
for large x and definition of r in (3.44c), it is obtained that 
(f)f = ^(0.0347r)-^/^[l+3.04xl0^^t]-^/2cos[0.068(l+3.04xl0^^t)-1.4997r]. 



(3.45a) 



So, 



(5(/)(t,x) = E(0.0347r)-^/2[l + 3.04 x 10^H]-^/^x 

cos[2.07 X lO^H - lA777r]e-'^-'' (3.456) 

It shows that S(f){t, x) oscillates about = 0c with decaying ampli- 
tude. Decay of amplitude of a classical field shows loss of energy from 
the system [42]. It means that we have two phenomena (i) production 
of elementary particles and anti-articles of these during inflation and 
(ii) at the end radiation of energy at the end of inflation. The radiated 
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energy from quintessence field heats up the universe and created par- 
ticles set in thermal equilibrium quickly if these particles have sufficient 
interaction. 

Phcise 2. Radiation-dominated phase 

It is found above that , during phase 1, the universe inflates for a 
very small fraction of a second and eqs.(3.12a) and (3.12c) show that 
Pde falls by 3 orders upto the end of inflation. The density for released 
dark energy is obtained as 



Pde(released) = Pde(O) - pM = 2.175 X 10^^GeV\ (3.46) 



using eqs.(3.12a) and (3.12c). 

It means that due to decaying quintessence field 0, DE emits as 
radiation and elementary particles having thermal equilibrium with 
this radiation. It is parallel to the idea introduced by Bronstein in 1933 
[6]. Here, it is proposed that this radiation causes cosmic microwave 
background radiation. So, 



from eqs.(3.46) and (3.47). 

We find another radiation term in Friedmann equation (2.19a) emerg- 
ing spontaneously from curvature along with DE term. 

The emitted radiation, due to decay of 0, sets dark radiation also in 
thermal equilibrium. 

Moreover, from (3.46) and (3.47), it is obtained that 




(3.47) 



Initial temperature for this radiation is obtained as 



Te = 1.35 X lO'^GeV 



(3.48) 



Prd(^c) > Pde{tc)- 
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Using the equation of state for emitted radiation = Prd/3 (with 
Prd and prd being the pressure and energy density respectively) conser- 
vation equation yields 

Prd = 5xl0-Vcr(-)', (3.49) 

using p°d = 5 X lO-^cr where = 2.33 x 0.68 x IQ-'^^GeV . Here Oq 
is the present value of the scale factor. 

Emitted radiation and produced particles (due to decay of quintes- 
sence DE) also contribute to the energy density in FE (2.19a). As a 
result, (2.19a) is modified as 



8nG 



Pde + 5X 10-Vcr(- 

\ a 



Pdrd + Pdm 



(3.50a) 



(3.506) 



where dark radiation density is given by 

8nG _ 0.05A 

—5— Pdrd — ^ 

3 a* 
being obtained from (2.22b) and (3.9b) for Wde- 

Pdm decreases as with the scale factor. So, in the early universe 
when a{t) is very small radiation terms dominate over matter. More- 
over, (3.49) shows dominance of radiation over DE term when t > tc- 
Thus (3.50a) is approximated as 
SttG 



4 6.27 X 10^^^ 
+ a* 



(3.50c) 



As dark radiation also has thermal equilibrium with radiation emitted 

due to dark energy decay , (3.50c) reduces to 



8nG 



5 X 10-V°r 



ao\4 6.27 xlO^^^ 
+ 



which integrates to 



a{t) ^a^[l + 2f3{t-tc)] 



5xl0-i7o^p) , 
(3.51) 



(3.52) 
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with 

In terms of eq.(3.51) can be re-written as 
So, 

.1 day _ 5 .W-y ^a,y ^ ^ ^ io--(5^)""° (3.53) 

KadcpJ (1 + WdcjPdc ^ ^ \aJ 

using eq. (2.27a) and wjc = —29/30. 
Eq.(3.53) is integrated to 



r /r;r,N 39/20 n :<:u/3y 

a{t) = ae 1 + 2.56 x lO'^o ( J (</> - (/.c) . (3.54) 
Connecting eqs. (2.28b) and (3.54), it is obtained that 

V{<t.) = 4.3 X ^)'''° \l + 2.56 X 10-°f^)''''% - 0. 



20/39 



-2/39 
(3.55) 



3.3. Phase 3. Matter-dominated phase 



Like radiation, we have two types of matter (i) dark matter given 
by Pdm in (2.19a), which is non-baryonic and (ii) baryonic matter 
formed by elementary particles, produced during phase I (inflation) 
through various processes of standard cosmology such as nucleosyn- 
thesis, baryosynthesis and recombination of hydrogen, which will be 
discussed in some other paper. Matter, so produced, contributes the 
density ~ 0.04p°j. to the present matter density p'^ ~ 0.27p°j. as p'^^ ~ 
0.23p0,. 

Matter dominates over radiation, when p^, > Pr, which implies that 

0.27> 5 X 10-^(^^) (3.56) 
using p^ and Pr from eqs. (2.23b) and (3.49) respectively. 
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Defining red-sliift as 

1 + ^ = ^, (3.57) 
a 

the inequality (3.27) yields that matter dominates over radiation at 

z < 5400. 



It is consistent with WMAP [29] , which gives that matter decouples 
from radiation at time = 386kyr = 1.85 x 10^''GeV~^ and red-shift 

= — - 1 = 1089. (3.58) 

So, we take the red-shift Zd — 1089, when matter dominates over 
radiation. In this case, on the effective Priedmann equation (2.24) is 
approximated as 

including the contribution of created matter. (3.59) integrates to 

a{t) = ad[l + 27(t - td^/^ (3.60) 



with 

7 



yM27Ho(—y^^ = 2.96 x 10-=^^GeV 

using eq.(3.58). 

In terms of 0, eq.(3.59) can be re- written as 

■ 1 da \ 2 , 3 



. d(f)/ 
So, 

\aa(p/ (1 -|- Wde)Pde ^ o / \a/ 
using eq. (2.27a) and Wde = —29/30. 
Eq.(3.32) is integrated to 

a(t) = ad[l + 3.54 x 10-^"((/) - (/)<i)]'°/'^ (3.62) 
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Connecting eqs. (2.28b) and (3.33), it is obtained that 



V{(f)) = 4.33 X 10-^^[1 + 5.5 X 10-^^{(f) - 0^)]-'/'^ (3.63) 



Further, formation of massive stars takes place during the matter- 
dominated phase of the universe. Large structure formation is not 
addressed in this work, but it is interesting to know the time ,tB, when 
the presence of massive stars is possible in this model. This time can be 
calculated from observations of Gamma-ray bursts (GRBs) taking place 
in sub-luminous star-forming host galaxies. GRBs show the existence 
of massive stars. According to previous observations, GRBs are found 
at high red-shift z ~ 1.6. But the recently observed GRB 050904 was 
found at red-shit zb — 6.llo!i2^ beating all previous records. In what 
follows, tB is evaluated using this recent observation. 

Connecting (3.57) and (3.60), zb is obtained as 



1 + = ^ = ^[1 + 4.44 X 10-''itB - U)]-'/^ (3.64) 



Using eq.(3.58) for ao/ad and zb given above, (3.64) yields Ib from 
822 Myr to 918 Myr after big-bang showing that massive stars were 
present by this time. 

3.4 Phase 4. Dominance of dark energy in the late universe 

Dominance of DE over matter (pde > Pm) is possible when 



\ a / 

using and Pr from eqs. (2.23a) and (2.23b) respectively. 

The inequahty (3.65) shows that dominance of DE, in the late uni- 
verse, is possible, when 



a(tB) ad 




(3.65) 



l + z< 1.409 



(3.66) 



a 



USmg Wde = - 



29/30 obtained above. 
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So, we can take red-shift for the transition from pde < Pm to pde > Pm 

as 

= 0.4 (3.67) 

safely. This is the red-shift at which DE dominates over matter. As 
cosmic acceleration is caused by dominance of DE. is the red-shift 
for transition from deceleration to acceleration. 

Interestingly, it falls in the range of red-shift 2; = 0.46 ± 0.13 for this 
transition, given by 16 Type supernova (SNe la) [5]. 

In this case, the Priedmann equation (2.24) reduces to 

(:)^=--"(7)"° 

which integrates to 

a{t) = a*4l + u{t - U^)f^, (3.69) 

with 

. = 0.04/^0 (^)''"(^)^'" = 6.88 X 10--GeV 

and a** = a{t^^). Here t** is the transition time. At the present time 
to = 13.7Gyr = 6.6 x lO^^GeV"^ we have 

1 + = = [1 + 6.88 X 10-^^(io - 

Connecting eqs.(3.67) and (3.69), it is evaluated that 

t„ = 4.13 X lO^^GeV"^ = S.SSGyr. (3.70) 

Eq.(3.69) shows accelerated expansion beginning at i**. So, 8.58Gyr 
is the time for transition from deceleration to acceleration in the late 
universe. 

In terms of 0, eq.(3.3) has the solution (for this case) as 

a = a,*e[^^^('^-'^**)l, (3.71) 
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where = 0(t„). 

Connecting eqs. (2.28b) and (3.71), it is obtained that 

V{(f)) = 2.22 X io37e[-V4^(0-</>..)]_ ^3 72) 

4. Universe driven by dark energy 

Eq.(3.12a) gives initial value of DE as 2.18 x lO^^GeV^ When it 
falls to its current value 2.18 x 10~''^GeV^, the ratio of Qq (present 
scale factor) and a(0) (the initial scale factor) is obtained as 10^^^° 
given by eq.(3.13b). It means that if universe is derived by dark energy 
only ao = 10^^^°a{0). 

In this case, cosmic dynamics is given by eq.(3.5) yielding the scale 
factor (3.6) and the expansion rate as 

^ - 7 ^ , / .1/20. (4-1) 

with Wde = —29/30 obtained above. Using eqs.(3.13a,b) in eq.(4.1) 

20 , , 

H (4.2) 

At t = to = 13.7G'yr = 6.6 x 10^^GeV"\ eq.(4.2) yields 

^0 = 3.05 X lO-^^GeV = 19.25//o- (4.3) 

It means that present Bubble's rate Hq ^ Hq = 19.25ifo in the 
universe driven by DE only. It happens due to highly accelerating 
expansion driven by DE having density (2.23a). The equation (4.3) 
shows that, in this situation, H is rescaled as 19.25if = H. With this 
rescaling, eq.(3.1b) looks like 
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and eq.(3.6) takes the form 

r 1 , / an \ 1/20 1 20 

a(t) = a(0)[l + -v^i/o(^) t\ , (4.5) 



yielding 



5.23 X 10^^^°. (4.6) 



a(0) a(0) 

This result is found for an universe, which is not being observed. We 
call it actual universe Now we come to our observable universe. 

4.1. Observable universe 

As discussed in the preceding section, universe passes through 4 
different phases ( power-law inflation, radiation-dominated, matter- 
dominated and DE dominated phase). We call it observable universe. 
According to analysis for observable universe, eq.(3.11) gives the ratio 

(^) =1.46x10^^. (4.7) 

At t — tc {tc being the time for end of inflation), eq.(3.49) yields 

15 X 10"^ ^^9/ao\^ 

Connecting eqs. (3.46)- (3.49) and (4.8), we obtain that 

(^) = 6.43 X 10^9. (4.9) 

\(lc/ obs. 

Prom eqs. (4. 7) and (4.9) 

9.39 X 10^^ (4.10) 



UfO)A 



l(0)/ obs 

4.2. Actual universe and Observable universe 

So, from eqs.(4.6)and (4.10), it is obtained that 

(^) ^h.2?,xW'»(^) (4.11) 

\a(Oj/ actual Va(0)/obs. 

The result (4.11) prompts us to think that actually, our observable 
universe resides in an extremely large universe driven by DE only. As 
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Pde(actual) 

= 370.56pde(obs.)- (4.12) 

It means that only ~ of DE is available in our observable universe. 

The reason for this scenario is given as follows. Due to the reverse 
gravity effect of DE, the actual universe speeds up very fast and accel- 
eration continues upto the present epoch giving the ratio (4.6), whereas 
the observable universe (driven by DE) inflates for a very short period 
in the beginning as well as decelerates driven by radiation and subse- 
quently by matter for the major portion of its age. Moreover, in the 
recent past, a transition from deceleration to acceleration takes place, 
which is consistent with 16 Type SNIa observations [5] Thus, the bi- 
furcation between the actual and observable universe begins in early 
times, when inflation of observable universe ends and its decelerated 
expansion, driven by radiation as well as dark matter followed by accel- 
eration in the recent past. As a result, for the observable universe, we 
have the ratio (4.10), which is extremely less the similar ratio for the 
actual universe given by (4.6). So, it is reasonable to conclude that the 
observable universe is concentrated in a very small part of the actual 
universe. 

5. Summary 

It is interesting to see that, in the present model, DE emerges sponta- 
neously due to manifestation of the Ricci scalar curvature i? as a phys- 
ical fleld too. The physical aspect of R is obtained adding a non-linear 
term to Einstein- Hilbert term. In [17, 18, 19, 20, 21, 22, 23, 24, 25], 
the non-hnear term of R stands for the lagrangian density of DE and 
it is added to explain late acceleration in the universe. In [20], Nojiri 



34 S.K.SRIVASTAVA 

and Odintsov have taken non-linear term ai?™ + with the motive 
to have inflation in the early universe and late acceleration, where R"^ 
yields acceleration in the early universe due to high curvature in the 
beginning of the universe and gives acceleration in the late uni- 
verse due to low curvature there. In contrast to this approach, here, 
only one non-linear term of R is taken in the action, which is not DE 
lagrangian. Here, DE is induced by linear and non-linear terms of R 
causing early and late acceleration. The resulting cosmological scenario 
is given as follows. In the beginning, our observable universe inflates 
for ~ 10~^^sec. DE, released during this period, decays to emission 
of CMB radiation and light elementary particles. This event heralds 
usual big-bang expansion with an end of inflation as emitted radiation 
dominates over DE by the time 10~^^sec. After this epoch, universe 
decelerates for a long time, driven by radiation and subsequently by 
matter. Around the age 8.58Gyr of the universe, a transition takes 
place and DE re-dominates over matter, giving a cosmic jerk. As a 
result, decelerating universe gets speeded up [?]. DE density, obtained 
here, falls from 2.58 x lO^^GeV^ to its current value 2.18 x 10-^^GeV^ 
Interestingly, it is found that the observable universe, having this sce- 
nario, expands within the actual universe, which is not observed. 

In non-gravitational field theoretic models of DE, a scalar field (pit) 
(quintessence, k-essence, tachyon and phantom ) are taken as DE source. 
In contrast to this approach, here, (f){t) giving DE originates from grav- 
itational sector. Thus, here, (f){t) has a geometrical origin, whereas 
origin of quintessence, k-essence and phantom is not known (except 
tachyon having stringy origin). The potential V{4>) ,for DE obtained 
here, is tabulated below in different situations. 
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Table no.l 
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Potential V{(j)) during different phases of tlie universe 



V{<j>) 



Dominance in Nature of 
the universe expansion 



2.56 X 10-20 (^^j (0_0^)]-2/39 

4.33 X 10-^^[1+ 
5.5 X 10-2=^(0 



7.12 X I062e-I'^-'^(0)1V¥^ 

/ \ 1/10 

4.3x10-4^^^) [1+ 




DE acceleration 



radiation 



deceleration 



DM 



deceleration 



2.22 X iQ^7e{-V^^^(<t>-<f>**)] 



DE 



acceleration 



In [47, 48], considering V{4>) ~ e~'^'^ (c being a constant) cosmic 
acceleration is derived in flat homogeneous model of the universe. Here 
the result is opposite. In the present paper, it is obtained that when 
universe accelerates , potential falls exponentially and it falls as power- 
law in case of deceleration. 



Boson case 

In the case of minimal coupling of $ to gravity and having interaction 
with (p as — (l/2)g(02$*$, we have the equation 



Appendix A 



□$(i,x) + +m|]$(i,x) = 



(Al) 



with m$ being mass of x). 
$ is decomposed as 



oo 



— oo 
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In (A2), are creation (annihilation) operators satisfying quan- 

tum conditions 

[ofe, ak'] = = [aj., 4,] 
[afe,4,] = {27rfSkk' (^3,4,5) 
The scalar product [41] is defined as 

($1, $2) = -i / ^/^d^x[^idt^ - {A6) 

J t=constant 

with dt denoting derivative with respect to time t . 

The 111— and out— states of $ are obtained at two extremes of the 
space-time where space-time is asymptotically Minkowskian. These are 
connected through Boglubov transformations 

This scalar product (A6) shows that = 0, ^aT) = 

1 = ($^«*, $^«*). So From (A7), it is obtained that 

obeying the condition 



\ak\' - ^ 1. (AW) 

Case of spin-1/2 Dirac fermion 

The equation of Dirac field ^' with mass is given as 

[ij^D^^ - mv,]* = (An) 

with i = y^,!?^ = + lr^^/i:)^.A7"7- Here ^ = = 
0,1,2,3), g^"" = KKt]"-^, {1^,1''} = 2^"'' and {7^,7^ = 277"^ with 
7'^(7") being Dirac matrices in curved(fiat) space-time. 
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Further ^ are decomposed in discrete modes k and spin (s/2) witli 
s = ±1 as 

oo 

s=±l k=—oo 
oo 

s=±l k=—oo 

with bk,s{blg) and d^k,s{d^-k s) being annihilation (creation) operators 
for positive (negative) energy particles respectively. Further, using 

*77fe,. = gkAvy''-%: (Alb) 

where 





/1\ 




/o\ 

















1 






















,Ui = 


1 





















vJ 



Here also, in- and out- states of ijj arc obtained at two extremes of 
the space-time where space-time is asymptotically Minkowskian. 
The Bogolubov transformations for i/j are given as 

id%-s)^ = C"M + ^-t-A. 
d^-t-s = ^sKj + Plsd'^k,-s- 
The Bogolubov coefficients a^^s and Pk,s satisfy the condition 

\ak,s\^ + \Pk,sf^l: {A16) 

where 

c^k,s = / d'x^\ls^Z'l iAl7) 

J t=constant 
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and 

J t=constant 
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